Kondo effect in the presence of itinerant-electron ferromagnetism 
studied with the numerical renormalization group method 
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The Kondo effect in quantum dots (QDs) - artificial magnetic impurities - attached to ferromag- 
netic leads is studied with the numerical renormalization group method. It is shown that the QD 
level is spin-split due to presence of ferromagnetic electrodes, leading to a suppression of the Kondo 
effect. We find that the Kondo effect can be restored by compensating this splitting with a magnetic 
field. Although the resulting Kondo resonance then has an unusual spin asymmetry with a reduced 
Kondo temperature, the ground state is still a locally-screened state, describable by Fermi liquid 
theory and a generalized Friedel sum rule, and transport at zero temperature is spin-independent. 

PACS numbers: 75.20.Hr, 72.15.Qm, 72.25.-b, 73.23.Hk 



Introduction. - The prediction and experimental 
observation of the Kondo effect in artificial magnetic im- 
purities - semiconducting quantum dots (QDs) 0, 0] - 
renewed interest in the Kondo effect and opened new 
opportunities of research. The successful observation of 
the Kondo effect in molecular QDs like carbon nanotubes 
0, |]| and single molecules attached to metallic elec- 
trodes opened the possibility to study the influence of 
many-body correlations in the leads (superconductivity 
or ferromagnetism) on the Kondo effect. Recently 
the question arose whether the Kondo effect in a QD at- 
tached to ferromagnetic leads can occur or not. Several 
authors have predicted |& |!J llOl LUj that the Kondo ef- 
fect should occur. However, it was shown recently 0] 
that the QD level will be spin-split due to the presence 
of ferromagnetic electrodes leading to a suppression of 
the Kondo effect, and that the Kondo resonance can 
be restored only by applying an external magnetic field. 
The analyses of Refs. |a IE LLfJ, llll ll2| were all based 
on approximate methods. 

In this Letter we resolve the controversy by adapting 
the numerical renormalization-group (NRG) technique 
[kH H3 | , one of the most accurate methods available to 
study strongly-correlated systems in the Kondo regime, 
to the case of a QD coupled to ferromagnetic leads with 
parallel magnetization directions. We find that in gen- 
eral the Kondo resonance is split, similar to the usual 
magnetic- field- induced splitting 0, 0|. However, we 
find that by appropriately tuning an external magnetic 
field, this splitting can be full y co mpensatedand and the 
Kondo effect can be restored jl7| (confirming Ref. H3)- 
We point out that precisely at this field the occupancy of 
the local level is the same for spin up and down, n-p = np , 
a fact that follows from the Friedel sum rule [23 • More- 
over, we show that the Kondo effect then has unusual 



properties such as a strong spin polarization of the Kondo 
resonance and, just as for ferromagnetic materials, for 
the density of states (DOS). Nevertheless, despite of the 
spin assymetries in the DOS of the QD and the leads, 
the symmetry in the occupancy n-\ — Tip implies that the 
system's ground state can be tuned to have a fully com- 
pensated local spin, in which case the QD conductance 
is found to be the same for each spin channel, G-p = Gp. 

The Model. - For ferromagnetic leads electron-electron 
interactions in the leads give rise to magnetic order and 
spin-dependent DOS p r -\{u>) =fi p r i(u>), r = L,R. Mag- 
netic order of typical band ferromagnets like Fe, Co, and 
Ni is mainly related to electron correlation effects in the 
relatively narrow 3d sub-bands, which only weakly hy- 
bridize with 4s and Ap bands [19(. We can assume that 
due to a strong spatial confinement of d electron orbitals, 
the contribution of electrons from d sub-bands to trans- 
port across the tunnel barrier can be neglected j^- fn 
such a situation the system can be modeled by nonin- 
ter acting [2l| s electrons, which are spin polarized due 
to the exchange interaction with uncompensated mag- 
netic moments of the completely localized d electrons. In 
mean-field approximation one can model this exchange 
interaction as an effective molecular field, which removes 
spin degeneracy in the system of nonintcracting conduct- 
ing electrons, leading to a spin-dependent DOS. The An- 
derson model (AM) for a QD with a single energy level 
£d, which is coupled to ferromagnetic leads, is given by 



H = } j t rka c\ ka c rka + e d ^ n° + u 

rk<7 u 

+ Y<( V r,kdUrM + Vr*kcl, k J a ) - 9 ^ B BS Z . (1) 
rka 

Here c r ka and d a (n„ — S a da) are the Fermi opera- 
tors for electrons with momentum k and spin er in the 
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leads (r = L/R), and in the QD, V r u is the tunnel- 
ing amplitude, S z = (n-f — hi)/2, and the last term 
is the Zeeman energy of the dot. In general, all in- 
formation about spin asymmetry in the leads can be 
modeled by the spin-dependent hybridization function 
T rrT (uj) = 7rJ] fc <5(^ - ^ka)V^ k = 7rp ro .(w)l/ r 2 , where 
Vr,k = const = V r and p riJ {io) is the spin-dependent DOS. 

In order to understand the Kondo physics of a QD 
attached to two identical ferromagnetic electrodes with 
parallel configurations, it suffices to study, instead of the 
above general model [Eq. Q], a simpler one, which cap- 
tures the same essential physics, namely the fact that 
r, r |(cj) 7^ Y r \\UJ) will generate an effective local magnetic 
field, which lifts the degeneracy of the local level (even 
for B = 0). A simple (but not unique) way of modeling 
this effect is to take the DOS in the leads to be constant 
and spin-independent, p rG (io) = p, the bandwidths to be 
equal £)f = D±, and lump all spin-dependence into the 
spin-dependent hybridization function, T ra (u)), which we 
take to be ^-independent, T r(T (u>) = T ra . By means of 
a unitary transformation the AM [Eq. HJ] can be 
mapped onto a model in which the correlated QD level 
couples only to one electron reservoir described by Fermi 
operators otka with strength T a — J2 r ^ra, 



H = 




(2) 



gfi B BS z 



Finally, we parametrize the spin-dependence of r CT in 
terms of a spin-polarization parameter P = (I\ — TA/T, 
by writing r T(i) ee \T(l ± P), where T = T T + T±. 

In the model of Eq. (J2J, we allowed for T| ^ but 
not for Z)f ^ _Dj, as would be appropriate for real fer- 
romagnets, whose spin-up and down bands always have 
a Stoner splitting AD = Df — D±, with typical values 
AD/H t < 20% (for Ni, Co, Fe). However, no essential 
physics is thereby lost, since the consequence of taking 
Z?T 7^ D± is the same as that of taking T| ^ Ti , namely 
to generate an effective local magnetic field [24j . 

The occurrence of the Kondo effect requires spin fluc- 
tuations in the dot as well as zero-energy spin-flip excita- 
tions in the leads. Indeed, a Stoner ferromagnet without 
full spin polarization —1 < P < 1 provides zero-energy 
Stoner excitations |25| , even in the presence of an exter- 
nal magnetic field. 

Method. - The model [Eq. J5J] can be treated by 
Wilson's NRG method. This method, with recent im- 
provements related to high-energy features and finite 
magnetic field 0, 0] , is a well-established method to 
study the Kondo impurity (QD) physics. It allows one 
to calculate the level occupation n a = (h a ) (a static 
property), the QD spin spectral function, Im xf( w ) — 
T {i®(t)([S z (t), S 2 (0)])}, where T denotes the Fourier 
transform, and the spin-resolved single-particle spec- 
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FIG. 1: Spin-dependent occupation of the dot level at (a) 
B = and (b) B = — O.ir, as a function of spin polarization 
P. (a) For B = 0, the condition n-\ — nj. only holds at 
P = 0. (b) For finite P it can be satisfied if a finite, fine- 
tuned magnetic field, B COTa p(P), is applied, whose dependence 
on P is shown in (c). As expected, it is approximately linear 
dlH. Here U = 0A2D, e d = -U/3, V = U/6, and T = 0. 
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FIG. 2: a) QD spectral function A(u) = A a (ui) for sev- 
eral values of spin polarization P; inset: expanded scale of 
the spectral function around eF. b) The spin-resolved spec- 
tral function for fixed P. For P — > — P, we have A a — > A-„. 
c) Compensation of the spin splitting by fine-tuning an exter- 
nal magnetic field. Parameters U, e<j, F and T as in Fig. 



tral density A a (u,T,B,P) = -±lm.g£ (w) for arbi- 
trary temperature T, magnetic field B and polarization 
P [where G^(ui) denotes a retarded Green function]. 
From this we can find the spin-resolved conductance 
G« = ^ II ^ F r ff /_ 00 00 ^(^)(-^) with f(co) de- 
noting the Fermi function and A = TLa/TRo- describing 
the asymmetry in the couplings. We choose A = 1 below. 

Generation and restoration of spin splitting. - In this 
Letter, we focus exclusively on the properties of the sys- 
tem at T = in the local moment - Kondo regime, where 
the total occupancy of the local level, n — J2a n <^ ~ 1- 
The occurence of charge fluctuations of course broadens 
and shifts the position of the QD levels (for both spin 
up and down), and hence changes their occupation. For 
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P 0, the charge fluctuations and hence level shifts and 
level occupations become spin- dependent, causing the QD 
level to split |l2| and the dot magnetization —n^ to be 
finite (Fig. [TJ. As a result, the Kondo resonance, too, is 
spin-split [2(3, and weakened (Fig. [2} , similarly to the 
effect of an applied mag netic field 15, 161. This means 
that Kondo correlations are reduced or even suppressed 
in the presence of ferromagnetic leads. However, for any 
fixed P, it is possible to compensate the splitting of the 
Kondo resonance [Fig. [3f c) and EKc)] by fine-tuning the 
magnetic field to an appropriate value, i? comp (P), defined 
as the field which maximizes the height of the Kondo res- 
onance. This field is found to depend linearly [2^ on P 
[Fig.QIc)] (as predicted in 0] for U — > 00). Remarkably, 
we also find (throughout the local moment regime) that 
at -Bcomp the local occupancies satisfy rif = nj, [Fig.^b)]. 
The fact that this occurs simultaneously with the disap- 
pearance of the Kondo resonance splitting suggests that 
the local spin is fully screened at P C omp- 

Spectral functions. - We henceforth fix the magnetic 
field at B = B comp (P). To learn more about the prop- 
erties of the corresponding ground state, we computed 
the spin spectral function Im{x^(u;)} for several values 
of P at T — (Fig. UJ). Its behavior is characteristic 
for the formation of a local Kondo singlet: as a func- 
tion of decreasing frequency, the spin spectral function 
shows a maximum at a frequency w max which we asso- 
ciate with the Kondo temperature [i.e. /csTk = hw max at 
B = Scomp (P)], and then decreases linearly with to, in- 
dicating the formation of the Fermi liquid state 0| ■ By 
determining Tk(P) (from u) max ) for different P- values, 
we find that Xk decreases with increasing P [Fig. Efb)]. 
For metals like Ni, Co, and Fe, where P = 0.24, 0.35, and 
0.40 respectively, the decrease of Xk is rather weak, so 
the Kondo effect should still be experimentally accessible. 
Remarkably, both Irn{xJ(w)} and the spectral function 
A a (u)) collapse rather well onto a universal curve if plot- 
ted in appropriate units [Figs.^a) and 03c)]. This indi- 
cates that an applied magnetic field £? C omp restores the 
universal behavior characteristic for the isotropic Kondo 
effect, in spite of the presence of spin-dependent cou- 
pling to the leads. Fig. QJa) shows that the amplitude 
of the Kondo resonance is strongly spin dependent for 
ferromagnetic leads, which is unusual and unique. The 
nature of this asymmetry is related to the asymmetry 
of the DOS in the leads, and its value is exactly pro- 
portional to ~ l/r CT . As a result the total conductance 
G a is not spin dependent [Fig. E{b)]. This indicates the 
robustness of the Kondo effect in this system: if the ex- 
ternal magnetic field has been tuned appropriately, it is 
able to compensate the presence of a spin asymmetry in 
the leads by creating a proper spin asymmetry in the dot 
spectral density, thereby conserving a fully compensated 
local spin and achieving perfect transmission. 

Friedel sum rule. - Further insights can be gained 
from the Friedel sum rule, an exact T = relation 



[29j that holds for arbitrary values of P and B [30j. 
The interacting Green's function can be expressed as 
13 G&tM = [w- e dCT + iT a - with spin- 

dependent edo- and T a ; the former due to Zeeman split- 
ting (e d(T = £d — 1/2 og^iBB) and latter due to the ferro- 
magnetic leads. Here S CT (w) denotes the spin-dependent 
self energy. Now, the Friedel sum rule |29j implies that 
at T = 0, the occupancy n a and spectral functions can 
be written as 



A CT (e F ) 



1 1 -1 
Oa/TT = tan 

Z 7T 

sin (7rn CT ) 



Edcr - e F + Sg (e F ) 

r„ 



, (3) 
(4) 



where T^(uj) = Re S CT (w), and </> CT (w) is the 
phase of G^(uj). Since Sf(e F ) ^ an 
equal spin occupation, n-j- — n±, is possible only 

Sf(ep))/r;, 



for (e dT - e F + £^(e F ))/r T = (e d | - e F 
which can be obtained only for an appropriate exter- 
nal magnetic field B = B comp . For the latter, in the 
local moment regime (n ~ 1) we have = ss 0.5, 
so that </>f = cf>i as 7r/2, which implies that the peaks 
of A-f and A± are aligned. Thus, the Friedel sum rule 
clarifies why the magnetic field B comp at which the split- 
ting of the Kondo resonance disappears, coincides with 
that for which n-j- = n^. For B = B comp , the spin- 
dependent amplitude A CT (e F ) of Eq. @ and the conduc- 
tance G„ ~ r <T A (T (e F ) agree well with the abovemen- 
tioned NRG results [Fig. 0(a), HV)]. 
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FIG. 3: a) The spin spectral function Im{x?(<^)} as a function 
of energy, b) Dependence of Tk on spin polarization P. The 
solid line shows the prediction from Ref. 12] for the functional 
form T K (P)/T K {0) = exp[C arctanh(P)/P], where the best 
fit is obtained for C = —5.98. Equation (6) from Ref. [l2T | with 
(PT+pJJo = (r/7r)[//[|e d |(t/ + e d )] would lead to C = -4.19. 
c) QD spectral function for several values of P. Parameters 
U, e d , T, and T are as in Fig. [I] B = B comp {P). 
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FIG. 4: (a) Spin resolved QD spectral function amplitude 
Au (lu — 0) at the Fermi level, and (b) the QD conductance 
Ga, as functions of the spin polarization P, for B = B conl p(P) 
and symmetric couplings (Tlct = r_R CT ), with U, td, T, and T 
as in Fig. implying n a ~ 0.5. The dashed line in (a) is 
1/(1 ± P) (Eq. © with n CT = 0.5). As expected, we find 
G CT = e 2 /h, with a numerical error less than 1%. 



In conclusion, applying the NRG technique, we showed 
that the Kondo effect in a QD attached to ferromag- 
netic leads is in general suppressed, because the latter 
induce a spin splitting of the QD level, which leads to 
an asymmetry in the occupancy rij ^ n±. Remarkably, 
the Kondo effect may nevertheless be restored by apply- 

tuned such that 



comp ? 



ing an external magnetic field B, 
the splitting of the Kondo resonance is compensated and 
the condition nj — n± is fulfilled. Although the Kondo 
resonance is strongly spin polarized, it then features a 
locally-screened state, a spin-independent conductance, 
and a Kondo temperature which decreases with increas- 
ing spin asymmetry. 
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Note added. - After submission of our paper, a 
preprint [M. S. Choi et al. cond-mat/0305107 studying 
a similar problem using the NRG technique appeared, 
with conclusions consistent with our's. 
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